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Time evolution of quantum systems is of interest in physics, in chemistry, and, more recently, in 
computer science. Quantum computers are suggested as one route to propagating quantum systems 
far more efficiently than ordinary numerical methods. In the past few years, researchers have 
actively been improving quantum simulation algorithms, especially those in second quantization. 
This work continues to advance the state-of-the-art by unifying several diverging approaches under 
a common framework. In particular, it highlights the similarities and differences of the first and 
second quantized algorithms which are usually presented in a distinct fashion. By combining aspects 
of the two approaches, this work moves towards an online second quantized algorithm operating 
within a single-Fock space. This paper also unifies a host of approaches to algorithmic quantum 
measurement by removing superficial differences. The aim of the effort is not only to give a high-level 
understanding of quantum simulation, but to move towards experimentally realizable algorithms 
with practical applications in chemistry and beyond. 


In the recent past, quantum information and quantum 
computing have had far reaching consequences for mod¬ 
ern quantum chemistry. Quantum information, on the 
one hand, has inspired powerful new numerical meth¬ 
ods. Entanglement based ansatzes collectively known as 
tensor network methods have become the primer method 
for obtaining ground states of one-dimensional and quasi- 
one-dimensional electronic system with promising appli¬ 
cations to condensed matter [1] and chemical system 
[2, 3]. On other hand, new quantum algorithms contain 
much promise for efficient paths to simulating quantum 
dynamics as well as learning about ground state behav¬ 
ior. Quantum simulation is the idea of using quantum 
computational devices for more efficient simulation [4]. 
Since the dynamics are simulated by a quantum system 
rather than calculated by a classical computer, quantum 
simulation often offers exponential advantage over clas¬ 
sical simulation [5] for the calculation of electronic ener¬ 
gies, reaction rates, correlation functions and molecular 
properties. 

Although a series of algorithmic improvements in sec¬ 
ond quantized quantum simulation have come out [6, 7] 
focused on improving the algorithm presented in [8], none 
of these articles have attempted to leverage the group 
theoretic insights presented in Ref. [9]. Here, the effort 
began in Ref. [9] is continued by attempting to construct 
a second quantized algorithm that acts only in a par¬ 
ticular Fock space. This has parallels with the locality 
arguments presented in Ref. [6] but directly builds upon 
the machinery used for first quantized simulations [10- 
12 ]- 

To arrive at an algorithm that operates within a partic¬ 
ular Fock space, various aspects of quantum simulations 
must be unified. Consistent with the goal of unifying 
approaches to quantum simulation, various algorithms 
for quantum measurement are compared. Despite the 
various articles on this topic spanning decades [10, 13- 


18], there are only superficial differences between these 
approaches. This is illustrated by showing equivalences 
among the various presentations of quantum measure¬ 
ment. 

The article begins with basic concepts from lattice 
quantum theory including the proper definitions of first 
and second quantization. In the following section, two 
canonical quantum algorithms are compared and con¬ 
trasted. This naturally leads to a new approach to second 
quantized simulation. In the penultimate section, vari¬ 
ous schemes for quantum measurement are united into 
a single schema. A summary and outlook towards the 
future ends the article. 


LATTICE QUANTUM THEORY 

Two ways to approach lattice systems are with first 
or second quantization. The underlying lattice sites are 
completely arbitrary and could be, for example, Slater 
type orbitals, Gaussian type orbitals, quadrature points 
from classical polynomials, or plane waves. For the en¬ 
suing discussion, the lattice selection is done by choosing 
the eigenbases of one-body operators. For the remainder 
of the article, N is the number of particles and M is the 
number of sites. 

In first quantization, one deals with a M N di¬ 
mensional complex tensor with N indices denoted as 
’fy{x\X 2 ...Xjv) € C with Xi £ This 

must be a completely antisymmetric tensor meaning 
\R(- • • i- ■ ■ j ■ ■ •) = — • • j ■ ■ ■ i ■ ■ •). Here capital let¬ 

ters are used to indicate collective indices i.e. X = 

(x\ • ■ ■ xjy). The wave function in bra-ket notion is: 

MM M M n 

A>(x-iX2-..x n ) \xi...x n ) = X *(*)!*> 

Xi X2 XN X 
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where vectors \x\X 2 ---Xn) are the many-particlc basis 
vectors of the M N dimension space. 

In second quantization, one has a listing of the unique 
elements of the first quantized tensor, that is a vector D 
with complex elements. The indexing of vector D is 
given by ordered 7V-tuples. One writes the wave function 
in second quantization as 

(w) 

l*>= E D K a^ K2 ...\n) = J2^ K k^\n) 

k 1 <k 2 <...<k n k 

( 2 ) 

The operators {di,aj : 1 < i,j < M} satisfy ebay = 
— a,jdi and aja; = 1 y — <baj where ljy = dtj whenever 
the underlying orbital basis is orthonormal. The vector 
|f2) is an arbitrary fixed vector in an 2 M dimensional 
space called the vacuum state and ]T] indicates that the 
sum is over ordered tuples. 

In lattice systems, first and second quantization are 
connected by local-basis transforms of the form b x = 
^ cijCjx • The anti-commutation relations between the 
two bases is given by [b y ,d]] + = C, y . To illustrate the 
connection, first define X = i=1 bi and X m = J][^ TO h 


c Jlb i(x ijt> - 

(3) 

Cjibi(Xljl) Cj 1 b 2 

(X 2 jt) 

(4) 

N 



E(— l) m+1 C jlbm {Xmj}) 
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Cjibi ‘ ‘ ‘ C jibN 


(6) 

Uhvbi ‘ ‘ ‘ Uhv&jv 




The final expression is called the Pliicker embedding of 
matrix C into an antisymmetric space [19]. Writing Cj X 
as <f>j ( x ), one sees that this is the standard Slater determi¬ 
nants from quantum chemistry up to the normalization 
factor 1/y/m which accounts for local basis functions in¬ 
dependently normalized. 

Hamiltonians 

In non-relativistic electronic quantum mechanics, the 
Hamiltonian consists of kinetic energy, two-electron 
Coulomb interactions, and a scalar field representing the 
electronic environment e.g. the nuclear charges. 

In a lattice basis corresponding to eigenfunctions of 
the potential operator, the potential energy operator 
is V = V r a\d r and the kinetic energy operator is 
T = '}2pqT pq a) p d q . Here, it is needed that the M x M 
matrix T is a fixed but arbitrary real symmetric matrix. 
The kinetic energy operator is complex if time reversal 


symmetry is broken. Define momentum space as the ba¬ 
sis where T is diagonal and let {b p } be the corresponding 
anti-symmetric operators. 

In the continuum limit, the change of basis between 
cti and bj is via the Fourier transform. However, in a 
discrete basis, this may not be the case due to various 
methods for approximating the kinetic energy operator. 
Still, the position and momentum spaces are related by 
some local unitary transform. 

An electron-electron interaction is written as W = 
■^Y^Wpqrsdjjdjjdrds. For the Coulomb interaction, this 
operator is diagonal in the position basis i.e. W = 
4 Wpqd^a^qdqdp. Moreover, the integral W pq depends 
only on the distance between sites p and q. Note that all 
of the Hamiltonian terms are permutationally invariant. 


UNIFIED QUANTUM SIMULATIONS 

In this section, we will unify the first and second quan¬ 
tized algorithms in a common framework and then create 
an partially online second quantized algorithm. By on¬ 
line, it is meant that the computation is done coherently 
in the quantum hardware. Before turning to the main 
thrust of the paper, it is worth mentioning algorithms 
simulating spare Hamiltonians following Refs. [20, 21] 
have recently been exploited in the context of quantum 
chemical simulations by Ref. [22]. The ideas presented 
below may also be applicable in this setting. 

In first quantized quantum simulations [10-12], the sys¬ 
tem is evolved under the operators W + V and T in the 
position and momentum basis, respectively. The quan¬ 
tum Fourier transform is used to efficiently convert be¬ 
tween position space and momentum space. The algo¬ 
rithms requires N\ogM qubits to store the wave func¬ 
tion and evolves the system only in a single Fock space. 
Antisymmetry is only imposed on the initial state which 
is stored in a binary representation whereby m qubits 
can store M = 2 m sites. In this work, this is the A1 
algorithm. 

In second quantized quantum simulations, the sys¬ 
tem is evolved under the Hamiltonian H = ^{T pq + 
V P q)aj,a q + ^'^2W pqrs dj ) c^ q d r d s . In the standard algo¬ 
rithms [7, 8], the integrals {V pq , T pq , W pqrs } are pre¬ 
computed classically and the wave function is stored in 
unary representation where M qubits are needed for M 
sites. The exponentially large Fock space (with a to¬ 
tal of J2k (10 = 2 M states) is only used for simulations 
that require changes in particle number such as super¬ 
conducting interactions approximated by cqaj + aj aj, 
grand canonical simulations with Hamiltonian terms 
b i 'l2k a k a ki or electron affinity/ionization calculations. 
When the particle number is fixed to N electrons, the 
evolution occurs in a subspace with only (^) states. A2 
is used to refer to this algorithm. 
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Online second quantized simulation algorithm 


The desirability of simulating within a single Fock 
space was highlighted in Ref. [9]. The first steps to¬ 
wards a second quantized algorithm achieving that ulti¬ 
mate goal is presented. The insights stems from contrast 
that scaling of algorithm A2 does not depend on the 
number of electrons being simulated whereas A1 does. 
The goal of simulating in a single Fock will not be fully 
achieved here, but similar to Ref. [6], a scaling of 0(M 2 ) 
is also obtained, but for very different reasons. Ref. [6] 
attempts to exploit locality of the local basis in order to 
show that the number of non-trivial Coulomb integrals 
scale as 0(M 2 );whereas, here the techniques from A1 are 
not completely adapted to second quantization. For this 
reason, the kinetic energy operator needs to be imple¬ 
mented in all Fock spaces giving 0(M 2 ) contributions. 
The spatial terms will scale with the number of electrons 
rather than the number of basis functions. 

In Al, the algorithm computes all quantities on-the- 
fly using phase-kickback of classically computable quan¬ 
tities. The key ideas is to work in the shared eigenbasis of 
the scalar potential and the two-body interaction term. 
Computing the Hamiltonian terms must be done coher¬ 
ently using quantum version of classical algorithms such 
as addition or multiplication [11, 12], At a high-level, the 
two-electron term is computed as e.g. follows: 
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In Eq. (8), two points of clarification are useful: (1) the 
position is stored in binary representation and (2) a quan¬ 
tum version of a classical computation is effected. For the 
implementation of the Coulomb interaction potential, the 
algorithm is predicated on addition, multiplication, and a 
Newton-Raphson method for implementing square roots 
leading to 0(b 2 N 2 ) scaling with b binary bits of precision 
used for computation [11], In Eq. (9), the phase kickback 
method is used [12, 16]. 

The interesting point uncovered here is that this com¬ 
putation can be easily adapted for second quantized algo¬ 
rithms. Indeed, since the Hamiltonian is permutationally 
invariant, the previous method can be directly adapted 
as, 
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Note that here the positions are stored in unary repre¬ 
sentation where the third bit corresponds to the third 


position and that counting pairs of fermions is not dif¬ 
ficult even when the particle number increases beyond 
two. 

One of the key differences between Al and A2, is the 
exploitation of Fourier transformation to efficiently im¬ 
plement the kinetic energy operator in its eigenbasis. In 
Al, this change of basis is done via the quantum Fourier 
transform. In A2, there are two obstacles: First, the 
quantum Fourier transform is conducted in the binary 
representation rather than the unary representation used 
in A2. Second, the coefficients of the Fourier transform 
must be anti-symmetrized as in Eq. (6). These difficulties 
are probably amendable using conversions from unary to 
binary, quantum Fourier transform, and the reversible 
anti-symmetrization algorithm [23]. However, this naive 
approach seems to defeat many of the key advantages of 
the A2. Accordingly, an elegant solution to performing 
local change of basis in second quantization is highly de¬ 
sirable in order to exploit the simplicity of the kinetic 
energy in the momentum basis. 

Despite the inability to efficiently convert between lo¬ 
cal bases in second quantized simulations, this paper pre¬ 
scribes an approach where by only the kinetic energy is 
implemented as in previous works [7, 8] but the poten¬ 
tial and the two-body part of the Hamiltonian are com¬ 
puted using their shared eigenbasis with the algorithm 
adapted from Al as presented here. This online version 
of the second quantized algorithm requires only the pre- 
computation of the kinetic energy operator. The online 
portions in the position basis can may be sped up by 
quantum pre-computation [12]. The detailed cost anal¬ 
ysis of the computations required by the online version 
was previously done in [12]. The full comparison of the 
merits of the standard second quantized algorithm and 
the online version introduced here will appear in the near 
future. 


UNIFIED QUANTUM MEASUREMENTS 

Thus far, the focus has been on the algorithms for 
propagating the quantum system. Now, attention turns 
to the task of measuring properties of the system. Sev¬ 
eral schemes are considered and shown to be equivalent. 
In order of appearance, the methods are von Neumann’s 
measurement apparatus [10, 13], Ramsey spectroscopy 
[14], phase estimation algorithm [15], phase kickback [16], 
adiabatic quantum simulation [17], and the energy spec¬ 
trum algorithm [18]. All of these schemes are the same 
up to superficial differences. 

The earliest example of the quantum measurement 
scheme was attributed to von Neumann by Ref. [10]. 
There one has an interaction Hamiltonian H m t = A® p 
where A = ]T] A Q |a)(a| is the observable to be measured. 
With initial state |a) ® |a; = 0) = |a) (g) aJ|Q), the evolu- 
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tion for a time 2n/M yields 


e -2mAp/M | a )| 0 ) = \ a ) e -2ni\ a /Mp\ x — — | a )| x = Aq ) 

(13) 

To obtain this result, the definitions from the Section 
are used with p = ^ kb\bk and change of basis be¬ 
tween position and momentum basis defined by Ckn = 
exp (— 2ni kn/M) / \[M. 

In the quantum information community, the most well 
known example of quantum measurement is the one given 
by [15]: 


|0> 

|a) 


-—A — b n 


U 


|a) 


(14) 

where H is the Hadamard transform defined as the 
Fourier transform over M = 2 sites. 

Next, analysis shows that Kitaev’s scheme and the von 
Neumann measurements are essentially the same. Ki¬ 
taev’s scheme requires implementing a controlled U. In 
equations, this is |0)(0|®l + |l)(l|(g)[/. Since this matrix 
is block diagonal, its generator is as well. The generator 
of the identity is just the zero matrix and the genera¬ 
tor of U is the Hamiltonian whose spectrum is of inter¬ 
est, A. The generator of the controlled unitary is then 
Flint = 11) (11 (g> A. This form is the same as is found in 
schemes of [17] and [18]. 

To fully illustrate the point, let us focus on the von 
Neumann scheme described above. Consider the case 
that there are only two sites i.e. M = 2. The kinetic 
energy operator for a two site system is proportional to 
the Hadamard (Fourier) transform of a z 


0 1 
1 0 


= —{ l 1 
>/2 U "I 


1 W 1 1 

V2 Vi -1 
(15) 

Shifting and rescaling the kinetic energy such that p = 
1(1 — a x ), one can rewrite the coupling Hamiltonian as 


H = p® A = (H®1) 


l-a z 


A )(H«1). (16) 


This Hamiltonian is precisely the generator of the Kitaev 
circuit. 

It remains to show that Ramsey spectroscopy [14] is 
equivalent to the other schemes. Ramsey spectroscopy 
is used to measure an unknown energy splitting. The 
protocol begins with a ir/2 pulse, followed by free evo¬ 
lution for time T, and ends with a final tt/2 pulse. The 
population is transfered to a phase difference and after 
some time evolution, the change in phase due to energy 
differences can be converted back into population. 

To show equivalence, the idea of “phase kickback,” 
which is yet another description of the von Neumann 
scheme, introduced by Ref. [16[ is utilized. Because the 
state |a) is an eigenvector of A, the evolution under A 


does nothing but impart a phase on the controlled regis¬ 
ter. In our case, this is what sets the qubit frequency cuq. 
This phase is said to be “kicked” back onto the control 
qubit (hence the name). From a reductionist point of 
view, the wave function register is superfluous and one 
can write the Kitaev circuit with a phase gate on the 
probe register: 


H 


G(A a t) 


H 




(17) 

Here G{9) = |0)(0| + exp(—i0)|l)(l|. To see the equiva¬ 
lence with the other schemes, consider the generator of 
such a phase gate Hq- For Hq, one needs (0|1 Fg| 0) = 0 
and <1 |£T g |1> = Aq,. This is nearly identical to the Ram¬ 
sey spectroscopic scheme. The minor difference is ac¬ 
counted for by the fact that the tt/2 pulse is not ex¬ 
actly the same as the Hadamard transformation. This is 
merely selecting a different approximation to the kinetic 
energy operator, but the outcomes are the same with the 
probability of being the |1) state being F(1 + cos(A Q t)) 
in both schemes. 

It should be noted that the exact number of mea¬ 
surements needed to estimate the value of a parame¬ 
ter [24-27] depends on how past measurements are uti¬ 
lized for the current measurement, what choice is made 
for the momentum operator, and how the inference is 
done. Thus, while the form of the measurements does 
not change, the usage of information obtained from mea¬ 
surements can be optimized. 

Further, note that the eigenvalues are always measured 
relative to something. It is never an absolute energy and 
this is a simple consequence of Fourier theory. Stated 
differently, the Hamiltonian operators can be unbounded 
but the unitary group is compact. This compactness 
means that you cannot see the full spectrum and that 
eigenvalues must eventually ’’wrap around.” The aliasing 
of the eigenvalues arises because the dynamics are unaf¬ 
fected by a shift of energy. Hence, there is a gauge degree 
for freedom which can be arbitrarily selected e.g. as all 
charges infinitely separated in electrostatics or all atoms 
completely separated in thermo-chemistry. A through 
discussion of aliasing in the context of quantum phase 
estimation is can be found in Ref. [8]. 

From Hamiltonian complexity results [28, 29], the 
quantum measurement of an eigenvalue is not an easy 
task. This is because Hamiltonians are exponentially 
larger than the problem size, hence there are an expo¬ 
nential number of eigenstates. If one wishes to scan the 
“absorption” spectrum by changing the probe’s energy 
as in the scheme of Ref. [18], the input states must have 
polynomially large overlap with the eigenstates for peaks 
to occur at eigenvalues. Otherwise, it will take an expo¬ 
nential number of measurements to see enough tunneling 
events to determine the transition energy. 
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OUTLOOK 

This article focused on understanding and combin¬ 
ing various aspects of existing quantum simulation algo¬ 
rithms. It was also shown that, despite some information- 
based improvements, the von Neumann scheme for quan¬ 
tum measurement has been re-discovered many times 
over the half century since its introduction. The intro¬ 
duction of a partially online second quantized algorithm 
opens the door to a fully online second quantized quan¬ 
tum simulation but this will require an intelligent method 
for changing the local basis from position space (eigen- 
basis of V) to momentum space (eigenbasis of T). 

Many experimental groups worldwide are pushing for 
more precise control over larger realizations of quantum 
computers, but simulations that enable quantum compu¬ 
tation to produce results unobtainable by standard nu¬ 
merical methods have yet to appear. As a model chem¬ 
istry [30], quantum simulation is still stuck in the third 
of the five stages of development: (1) target accuracy, (2) 
formulation, (3) implementation, (4) validation, and (5) 
prediction. 

The current work continues to simplify the formulation 
in an effort to make the implementation more feasible for 
experimental technologies. But to move on to the fourth 
stage, a test set of molecular instances must be tested in 
functional quantum hardware. Until then, the ultimate 
goal of making predictions will remain far into the future. 

Taking an eye towards near-future schemes for quan¬ 
tum simulation, a paradigm shift from competition to 
collaboration is needed. Instead of asking how a quan¬ 
tum computer can beat a classical computer, we should 
be asking how a quantum sub-processor can enhance 
commonplace processors. The necessary developments 
to bring this to fruition are (1) applications where small 
quantum processors can outperform small CPUs, (2) ap¬ 
plications that require small, short quantum evolutions 
as part of a larger application (e.g. crucial quantum 
coherent subsystems) and (3) a quantum-classical data 
bus. These are the components that will make quantum 
computation a viable technology that will become part 
of commercial products. This shifts the emphasis from 
larger quantum computers to faster, more reliable ones 
and that should only require modest advances in experi¬ 
mental technology. 
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